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(Dated: May 21, 2001) 
We present measurements of the branching fractions and a search for CP-violating charge 
asymmetries in charmless hadronic decays of B mesons into two-body final states of kaons and 
pions. The results are based on a data sample of approximately 23 million BB pairs collected 
by the BABAR detector at the PEP-II asymmetric B Factory at SLAC. We find the following 
branching fractions: B(B° -> 7r + 7r~) = (4.1 ± 1.0±0.7) x 1(T 6 , B(B° -» K+ty') = (16.7±1.6± 
1.3)xl(T 6 , B(B+ ^K+tv°) = (10.8ttl ±1.0) xHT 6 , 23(5+ -> K°tv+) = (18.2±|;g ±2.0) x HT 6 , 
B(B° -» K°ty°) = (8.2±|;^ ± 1.2) x 10" 6 . We also report the 90% confidence level upper limits 
B(B° K+K-) < 2.5xl0~ 6 , B(B+ -> tt+tt ) < 9.6xl0~ 6 , and B(B+ -» K°K+) < 2.4xl0~ 6 . 
In addition, charge asymmetries have been measured and found to be consistent with zero, where 
the statistical precision is in the range of ±0.10 to ±0.18, depending on the decay mode. 

PACS numbers: 13.25.Hw, 13.25.-k, 14.40.Nd 



The study of B meson decays into charmless hadronic 
final states plays an important role in the understand- 
ing of CP violation. In the Standard Model, all CP- 
violating phenomena are a consequence of a single com- 
plex phase in the Cabibbo-Kobayashi-Maskawa (CKM) 
quark-mixing matrix 0. Recently, the Belle and BABAR 
collaborations published results |^] on measurements 
of CP- violating asymmetries in B decays into final states 
containing charmonium, leading to constraints on the an- 
gle f3 of the CKM Unitarity Triangle. Measurements of 
the rates and charge asymmetries for B decays into the 
charmless final states 7T7t and Kit can be used to con- 
strain the angles a and 7 Q of the Unitarity Triangle. 

In this Letter we present new measurements of the 
branching fractions for B meson decays to the charm- 
less hadronic final states 7r + 7r _ , K + ir~, K + tt°, K°tt + 
and K°ir° ||. In addition, we search for charge asym- 
metries in the modes B° — > K + ir~, B + — > K + ir° and 
B + — > K°tt + . Previous measurements |], [tJ of these 
decays were reported by the CLEO Collaboration. 

The data sample used in these analyses was collected 
with the BABAR detector [g at the PEP-II e+e~ collider 
H at the Stanford Linear Accelerator Center. It cor- 
responds to an integrated luminosity of 20.6 fb _1 taken 
on the T(4S) resonance ( "on-resonance" ) and 2.61 fb _1 
taken at a center-of-mass (CM) energy 40 MeV below the 
T(4S) resonance ( "off- resonance" ) , which are used for 
continuum background studies. The on-resonance sam- 
ple corresponds to (22.57 ± 0.36) x 10 6 BB pairs. The 
collider is operated with asymmetric beam energies, pro- 
ducing a boost (/?7 = 0.56) of the T(45) along the colli- 
sion axis (z). The boost increases the momentum range 
of two-body B decay products from a narrow distribution 
centered near 2.6 GeV/c to a broad distribution extending 
from 1.7 to 4.3GeV/c. 

The BABAR detector is a spectrometer of charged and 
neutral particles and is described in detail in Ref. ||. 
Charged particle (track) momenta are measured in a 



tracking system consisting of a 5-layer, double-sided, sil- 
icon vertex detector and a 40-layer drift chamber (DCH) 
filled with a gas mixture of helium (80%) and isobutane 
(20%), both operating within a 1.5 T superconducting 
solenoidal magnet. Photons are detected in an electro- 
magnetic calorimeter (EMC) consisting of 6580 CsI(Tl) 
crystals. Charged hadron identification is based on the 
Cherenkov angle 9 C measured by a unique, internally re- 
flecting Cherenkov ring imaging detector (DIRC). 

Hadronic events are selected based on track multiplic- 
ity and event topology. Backgrounds from non-hadronic 
events are reduced by requiring the ratio of Fox- Wolfram 
moments H2/H0 |l(]] to be less than 0.95 and the spheric- 
ity jllj of the event to be greater than 0.01. 

All tracks (except K® decay products) are required to 
have a polar angle within the tracking fiducial region 
0.41 < 6 < 2.54 rad and a Cherenkov measurement from 
the DIRC. The latter is satisfied by 91% of the tracks 
in the fiducial region. We require a minimum number of 
Cherenkov photons associated with each 9 C measurement 
in order to improve the resolution. The efficiency of this 
requirement is 97% per track. Tracks with a 9 C within 3<7 
of the expected value for a proton are rejected. Electrons 
are rejected based on specific ionization (dE/dx) in the 
DCH system, shower shape in the EMC, and the ratio of 
shower energy to track momentum. 

Candidate mesons are reconstructed from pairs of 
oppositely charged tracks that form a well-measured ver- 
tex and have an invariant mass within 3.5<7 of the nominal 
Kg mass jL2j. The measured proper decay time of the 
K^ candidate is required to exceed 5 times its error. 

Candidate ir° mesons are formed from pairs of photons 
with an invariant mass within 3cr of the nominal ir° mass. 
Photons are defined as showers in the EMC that have the 
expected lateral shape, are not matched to a track, and 
have a minimum energy of 30 MeV. The ir° candidates 
are then kinematically fitted with their mass constrained 
to the nominal ir° mass. 
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B meson candidates are reconstructed in four topolo- 
h+h'-, h+ir°, K° s h+ and K° s ir°, where the sym- 
h and h! refer to 7r or K. The kinematic con- 
straints provided by the T(4S) initial state and rela- 
tively precise knowledge of the beam energies are ex- 
ploited to efficiently identify B candidates. We define a 
beam-energy substituted mass toes = \J E^ — p 2 B , where 
Eh = (s/2 + pi-pB)/Ei, ^/s and Ei are the total energies 
of the e + e~ system in the CM and lab frames, respec- 
tively and pi and ps are the momentum vectors in the 
lab frame of the e + e~ system and the B candidate, re- 
spectively. To improve the resolution in modes contain- 
ing 7T° mesons, the B candidate is kinematically fitted 
with the energy constrained to the CM beam energy For 
all modes, the toes resolution is dominated by the beam 
energy spread and is approximately 2.5MeV/c 2 . Candi- 
dates are selected in the range 5.2 < toes < 5.3GeV/c 2 . 

Wc define an additional kinematic parameter AE as 
the difference between the energy of the B candidate 
and half the energy of the e + e~ system, computed in 
the CM system, where the pion mass is assumed for all 
charged decay products of the B. The AE distribution 
is peaked near zero for modes with no charged kaons 
and shifted on average — 45MeV (— 91 MeV) for modes 
with one (two) kaons, where the exact separation de- 
pends on the laboratory kaon momentum. The resolution 
on AE is mode dependent. For final states that contain 
no 7T° mesons the resolution is about 26 MeV. For modes 
with 7T° mesons the resolution is about 42 MeV and is 
asymmetric due to underestimation of the ir° energy in 
the EMC. Candidates are accepted in the following AE 
ranges (given in GeV): [-0.15,0.15] {h+h'~), [-0.2,0.15] 
(h+n"), [-0.115,0.075] (K°h+) and [-0.2,0.2] {K° s n°). 

Detailed Monte Carlo simulation, off-resonance data, 
and events in on-resonance Toes and AE sideband regions 
are used to study backgrounds. The contribution due to 
other £>-meson decays, both from b — > c and charmless 
decays, is found to be negligible. The largest source of 
background is from random combinations of tracks and 
neutrals produced in the e + e~ — ► qq continuum (where 
q = u, d, s or c). In the CM frame this background 
typically exhibits a two-jet structure that can produce 
two high momentum, nearly back-to-back particles, in 
contrast to the spherically symmetric nature of the low 
momentum T(4S) — ► BB events. 

We exploit this topology difference by making use of 
two event-shape quantities. The first variable is the angle 
S [[llj between the sphericity axes of the B candidate 
and of the remaining tracks and photons in the event. 
The distribution of | cos 9 S | in the CM frame is strongly 
peaked near 1 for continuum events and is approximately 
uniform for BB events. We require | cos6> s | < 0.9, which 
rejects 66% of the background that remains at this stage 
of the analysis. 

The second quantity is a Fisher discriminant T con- 
structed from the scalar sum of the CM momenta of 



all tracks and photons (excluding the B candidate de- 
cay products) flowing into nine concentric cones cen- 
tered on the thrust axis of the B candidate. Each 
cone subtends an angle of 10° and is folded to com- 
bine the forward and backward intervals. Monte Carlo 
samples are used to obtain the values of the coefficients, 
which are chosen to maximize the statistical separation 
between signal and background events. The distribu- 
tions of T for Monte Carlo simulated B° — + h + h'~ de- 
cays and background events in the toes sideband region 
5.20 < toes < 5.27GeV/c 2 are displayed in Fig. [11(a). 
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FIG. 1: (a) The distributions of the Fisher discriminant 
for Monte Carlo simulated B° — > h + h'~ decays (histogram) 
and background events (points) in the wies sideband region 
5.20 < mBS < 5.27 GeV/ c 2 ; (b) the K-n separation, in units 
of standard deviations, as a function of momentum, derived 
from the Cherenkov angle measurements of kaon and pion 
tracks in a D* + — > D°ir + control sample, as described in the 
text. 



The final reconstruction efficiencies range from 31% to 
45%, depending on the mode. The detection efficiencies, 
which include the branching fractions of K° — > — > 
7r + 7r~ and 7T° — > 77 [Q, are listed in Table Q. 

Signal yields are determined from an unbinned maxi- 
mum likelihood fit that uses toesj AE, !F, and 6 C (where 
applicable). Separate fits are performed for each of the 
four topologies, where the likelihood for a given candidate 
j is obtained by summing the product of event yield u, 
and probability Vi over all the possible signal and back- 
ground hypotheses i. The rij are determined by maxi- 
mizing the extended likelihood function C: 



M 



N 



C = exp 



& n 



(1) 



The probabilities Vi(xj)oii) are evaluated as the prod- 
uct of probability density functions (PDFs) for each of 
the independent variables Xj, given the set of parame- 
ters d?j. Monte Carlo simulation is used to validate the 
assumption that the fit variables are uncorrelated. The 
exponential factor in the likelihood accounts for Poisson 
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TABLE I: Summary of results for detection efficiencies (e), fitted signal yields {Ns), statistical significances (5), measured 
branching fractions (B), and charge asymmetries. The efficiencies include the branching fractions for K° — > — > tt + tt~ and 
ty° — > 77. Equal branching fractions for T(4S) — > B°B° and B + B~ are assumed. The 90% confidence level (C.L.) intervals for 
the charge asymmetries include the systematic uncertainties, which have been added in quadrature with the statistical errors. 



Mode 


£ <%) 


1 v 


S (a) 




A 


A 90% C L 


K+TT- 

K+K~ 


45 
45 
43 


41 ± 10 ± 7 
169 ± 17 ± 13 
8.2±H ±3.5 


4.7 
15.8 
1.3 


4.1 ± 1.0 ± 0.7 
16.7 ±1.6 ±1.3 
< 2.5 (90% C.L.) 


-0.19 ±0.10 ±0.03 


[-0.35, -0.03] 


n+n 
K+n° 


32 
31 


37 ± 14 ±6 
75 ± 14 ± 7 


3.4 
8.0 


< 9.6 (90% C.L.) 
10.8±?;£ ± 1.0 


0.00 ±0.18 ±0.04 


[-0.30, ±0.30] 


K°n+ 
K°K+ 


14 
14 


59llo ± 6 
-4.1±|;| ± 2.3 


9.8 


18.2l^ ± 2.0 
< 2.4 (90% C.L.) 


-0.21 ±0.18 ±0.03 


[-0.51, ±0.09] 


K°7T° 


10 


17.9±g J ± 1.9 


4.5 


8.2tf;£ ± 1.2 







fluctuations in the total number of observed events N. 
For the K+tt^, ti^tt , K ± 7r°, J ftT°7r ± , and if°K ± terms, 
the yields are rewritten in terms of the sum nf + nj and 
the asymmetry A = (nj—nf)/ (nj+nj), where rif (nj) is 
the fitted number of events in the mode B — > / (B — *■ /). 
The numbers of events, N, entering the maximum like- 
lihood fit for each topology are 16032 (h+h'-), 16452 
(h+7r°), 3623 (K° s h+), and 1503 {K° s Tr°). 

The parameters for background toes and AE PDFs 
are determined from events in on-resonance AE side- 
band regions. The signal toes and AE PDF parameters 
are determined from fully reconstructed B + — ► D°tt + 
and B + — > D°p + (p + — > 7r + 7r°) decays. Events in on- 
resonance toes sideband regions and Monte Carlo simu- 
lated signal decays are used to parameterize the Fisher 
discriminant PDFs for background and signal, respec- 
tively (see Fig. 0(a)). Alternative parameterizations ob- 
tained from off-resonance data and Monte Carlo simu- 
lation are used as cross-checks and for determination of 
systematic uncertainties. The 8 C PDFs are derived from 
kaon and pion tracks in the momentum range of interest 
from approximately 42 000 D*+ -> D°ir+ (D° -> K~Tr + ) 
decays. This control sample is used to parameterize 
the Be resolution ae a as a function of track polar angle. 
The resulting K-ir separation, defined as \8f — d^\/ag c: 
where Of (0^) is the expected Cherenkov angle for a kaon 
(pion), is shown as a function of momentum in Fig. |l|(b). 

The results of the fit are summarized in Table |, where 
the statistical error for each mode corresponds to a 68% 
confidence interval and is given by the change in signal 
yield rii that corresponds to a — 21n£ increase of one 
unit. Signal significance is defined as the square root of 
the change in — 2 ln£ with the corresponding signal yield 
fixed to zero. For the three modes that have statistical 
significance less than 4<r we report Bayesian 90% confi- 
dence level upper limits. In addition, for the purpose of 
combining with measurements from other experiments, 



we report the branching fractions corresponding to the 
fitted signal yields: B{B+ -> tt+tt ) = (5.1±^ ± 0.8) x 
10- 6 , B{B° K+K-) = (0.85±g;|^ ± 0.37) x 10~ 6 and 
B(B+ -> K°K+) = (-1.3^ ± 0.7) x lO" 6 . The upper 
limit on the signal yield for mode i is given by the value 

of n° for which J™* £ max dm/ J °° £ max drii = 0.90, where 
jC m ax is the likelihood as a function of rii , maximized with 
respect to the remaining fit parameters. Branching frac- 
tion upper limits are calculated by increasing the signal 
yield upper limit and reducing the efficiency by their re- 
spective systematic errors. 

Figure || shows the distributions in toes and AE for 
events passing the selection criteria, as well as require- 
ments on likelihood ratios, which are used to increase 
the relative fraction of signal events of a given type. 
These likelihood ratios are defined for a given topology 
as K sig = J2s n s'Ps/ J2i n iPi and T^k = n k V k / XL n s V s , 
where ^ s denotes the sum over the probabilities for sig- 
nal hypotheses only, J2i denotes the sum over all the 
probabilities (signal and background), and Vk denotes 
the probability for signal hypothesis k. These probabil- 
ities are constructed from all the PDFs except that de- 
scribing the displayed variable. The likelihood fit projec- 
tions, scaled by the relative efficiencies for the likelihood 
ratio requirements, are overlaid on each distribution. 

Systematic uncertainties arise from: imperfect knowl- 
edge of the PDF shapes, uncertainties in the detection 
efficiencies, and potential charge bias in track reconstruc- 
tion and particle identification. Uncertainties in the PDF 
shapes affect both branching fraction and charge asym- 
metry measurements. 

For most of the branching fraction measurements, the 
PDF shapes contribute the largest systematic error. The 
exception is the B + — > K + ir° mode, where the largest 
systematic error is due to the 5% uncertainty in the ir° 
reconstruction efficiency. PDF systematic errors are es- 
timated either by varying the PDF parameters within 
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FIG. 2: The m,ES and AE distributions for the various 
modes, using likelihood ratio requirements described in the 
text. The solid curves represent the fit predictions for both 
signal and background; the dashed curve represents the given 
signal mode only and the dotted curve represents other modes 
of the same topology. 



tics samples of charged tracks in multihadron events. 
These studies show differences in reconstruction efficien- 
cies for positively and negatively charged tracks of less 
than 0.005. We assign an overall systematic uncertainty 
of ±0.01 on A for possible charge-correlated biases in 
track reconstruction and particle identification. All mea- 
sured background asymmetries are consistent with zero 
with statistical uncertainties less than 0.03. The fitted 
signal yields and asymmetries for off-resonance data and 
on-resonance AE sidebands are also consistent with zero. 

The overall systematic errors on the branching frac- 
tions and charge asymmetry measurements are computed 
by adding in quadrature the PDF systematic uncertain- 
ties and the systematic uncertainties on the efficiencies 
or due to possible charge biases, respectively. 

In summary, we have measured branching fractions for 
the rare charmless decays B° — ► 7r + 7r~, B° — ► K + ir~, 
B+ -> K+TT , B+ -> K tt+, and B° -> K°ir°, and 
set upper limits on B° — > K + K~, B + — > 7r + 7r°, and 
B + — > K K + ■ We find no evidence for direct CP viola- 
tion in the observed decays and set 90% C.L. intervals. 
These measurements are in good agreement with previ- 
ous results H 0] . 
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1 a of their measured uncertainties or by substituting al- 
ternative PDFs from independent control samples. The 
systematic errors in the signal yields due to PDF uncer- 
tainties depend on decay mode as shown in Table Q. 

The D* + control sample of kaon and pion tracks is used 
to estimate systematic uncertainties in the asymmetries 
arising from possible charge biases in the 6 C quality re- 
quirements, as well as from differences in 9 C reconstruc- 
tion for different charge species. From these studies we 
conservatively assign a systematic uncertainty of ±0.01 
on A for all the modes. Charge biases in the detector 
and track reconstruction chain are studied in high statis- 



* Also with Universita di Perugia, Perugia, Italy, 
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